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Abstract 

Total reaction cross sections have been predicted for nucleons scattering from nuclei ranging in 
mass 6 to 238 and for projectile energies from just above noticeable giant resonance excitation to 
300 MeV. So also have been the mass variations of those cross sections at selected energies when 
they have been calculated using coordinate space optical potentials formed by full folding effective 
two-nucleon (NN) interactions with one body density matrix elements (OBDME) of the nuclear 
ground states. Good comparisons with data result when effective A^A^ interactions defined by 
medium modification of free NN t matrices are used. However, there is a simple three parameter 
functional form that reproduces the partial wave total reaction cross section values determined 
from those optical potential calculations; a functional form also maps the total scattering cross 
section partial wave elements. Adjusting the theoretical defined parameter values has enabled us 
to fit the actual measured data values from the scattering involving (15) nuclei spanning the mass 
range from ^Be to ^^^U and for proton energies from 10 to 300 MeV. Likewise total cross sections 
for neutron cross sections for neutron scattering from various nuclei can be equally well reproduced. 
Of import is that the three parameter values vary smoothly with mass and energy. 
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I. INTRODUCTION 



Total and total reaction cross sections from the scattering of nucleons by nuclei are 
required in a number of fields of study which range over problems in basic science as well 
as many of applied nature. Such as in radioactive waste management, in radiation therapy 
and radiation protection for patients and in many other fields. In medical radiotherapy 
absorbed dose distributions in the patient are needed and cannot be measured directly as 
we should not do nuclear experiment with humen body, they must be calculated. For all 
these purposes an extensive data bank is necessary. Since we do not have the experimental 
data for all necessary nuclei at all different energies it would be utilitarian if such scattering 
data were well approximated by a simple convenient function form with which predictions 
could be made for cases of energies and/or masses as yet to be measured. It has been 
shown 0, 0, 0, 0] that such forms may exist for proton total reaction cross sections. Herein 
we consider that concept further to reproduce the measured total reaction cross sections 
from protons at energies from 10 to 300 MeV and from 15 different nuclei ranging in mass 
between ^Be and ^38^ g^j^^ total cross sections from neutron scattering for energies to 600 
MeV and from nine nuclei ranging in mass between ^Li and ^^^U. These suffice to show that 
such forms will also be applicable in dealing with other stable nuclei since their neutron 
total cross sections vary so similarly with energy 

Total reaction cross sections for protons and total cross sections for neutrons from nuclei 
have been well reproduced by using optical potentials. In particular, the data (to 300 MeV) 
from the same nuclei we consider, compare quite well with predictions made using a gf-folding 
method to form nonlocal optical potentials p , though there are some notable discrepancies. 
Alternatively, in a recent study Koning and Delaroche [s*] gave a detailed specification of 
phenomenological global optical model potentials determined by fits to quite a vast amount 
of data, and in particular to the neutron total scattering cross sections we consider herein. 



II. FORMALISM 

The total and total reaction cross sections for nucleons scattering from nuclei can be 
expressed in terms of partial wave scattering (5*) matrices specified at energies E o: k"^, hj 

^ S^ik) = e'^'^^'^ = r/±(A;)e^*K^('=)] , (1) 

where Sf^{k) are the (complex) scattering phase shifts and ri^{k) are the moduli of the S 
matrices. The superscript designates j = / ± 1/2. In terms of these quantities, the elastic, 
reaction (absorption), and total cross sections respectively are given by 
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= {l + l){l-vt{k)cos{2^[6t{k)\)] + l{l-r^;{k)cos[2^[6i{k)\)] .(4) 
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Therein the a\ are defined as partial cross sections of the total elastic, total reaction, and 
total scattering itself. For proton scattering, because Coulomb amplitudes diverge at zero 
degree scattering, only total reaction cross sections are measured. Nonetheless study of such 
data 0,0] established that partial total reaction cross sections a\^\E) may be described by 
the simple function form 



aj^) {E) = {21 + 1) [1 + J + e (2/o + 1) [l + J , (5) 

with the tabulated values of lo{E, A), a{E, A), and e{E, A) all varying smoothly with energy 
and mass. Those studies were initiated with the partial reaction cross sections determined 
by using complex, non-local, energy-dependent, optical potentials generated from a (7- folding 
formalism While those (7- folding calculations did not always give excellent reproduction 
of the measured data (from ~ 20 to 300 MeV for which one may assume that the method 
of analysis is credible), they did show a pattern for the partial reaction cross sections that 
suggest the simple function form given in Eq. (0). With that form excellent reproduction of 
the proton total reaction cross sections for many targets and over a wide range of energies 
were found with parameter values that varied smoothly with energy and mass. 

Herein we establish that the partial reaction cross sections for scattering of protons and 
total cross sections for scattering of neutrons from nuclei can also be so expressed and we 
suggest forms, at least first average result forms, for the characteristic energy and mass 
variations of the three parameters involved. Fifteen nuclei, ^Be, ^^C, ^^O, ^^F, ^''Al, ^°Ca, 
63Cu, 9°Zr, "^Sn, i^oq^^ i59rj^^ isi^g^^ iqt^u^ aosp^ ^^^^ 238u^ are considered for proton 

case and nine nuclei, ^Li, ^^c, i^F, ^°Ca, ^^Y, ^^^W, ^^^Au, ^ospb and ^ssu, for which a 
large set of experimental data exist, are considered for neutron case. Also those nuclei span 
essentially the whole range of target mass. However, to set up an appropriate simple function 
form, initial partial total cross sections must be defined by some method that is physically 
reasonable. Thereafter the measured total cross-section values themselves can be used to 
tune details, and of the parameter Iq in particular. We chose to use results from (7-folding 
optical potential calculations to give those starting values. 
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III. RESULTS AND DISCUSSIONS 

The function form results we display in the following set of figures were obtained by 
starting with gf-folding model results, for proton case at energies of 10 to 100 MeV in steps 
of 10 MeV, then to 300 MeV in steps of 50 MeV, and for neutron case, at energies of 10 
to 100 MeV in steps of 10 MeV, then to 350 MeV in steps of 25 MeV, and thereafter in 
steps of 50 MeV to 600 MeV. The (yf-matrices used above pion threshold were those obtained 
from an optical potential correction to the BonnB force which, while approximating 
the effects of resonance terms such as virtual excitation of the A, may still be somewhat 
inadequate for use in nucleon-nucleus scattering above 300 MeV. Also relativistic effects 
in scattering, other than simply the use of relativistic kinematics in the distorted wave 
approximation (DWA) approach, are to be expected. Nonetheless the DWA results are used 
only to find a sensible starting set of the function form parameters Iq, a, and e from which to 
find ones that reproduce the measured total cross-section data. One must also note that the 
(7-folding potentials for most of the nuclei considered were formed using extremely simple 
model prescriptions of their ground states. A previous study revealed that with good 
spectroscopy the ^f-folding approach gives much better results in comparison with data than 
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that approach did when simple packed shell prescriptions for the structure of targets were 
used. That was also the case when scattering from exotic, so-called nucleon halo, nuclei 
were studied 0, 0| • 

Although on using Eqs. (jH)) and (0) to match values of (theoretically) calculated total 
reaction cross sections led [2| to the three parameters lo{E, A), a{E, A), and e{E, A) having 
smooth variations with energy, there are discrepancies between those predictions and the 
actual measured data. Herein we improve the method of selection of those parameter values 
to produce more accurate reaction cross section values, while keeping as smooth a variation 
with energy of those parameters as possible. Specifically, in Eq. (jSj), we have set the e as a 
constant (—1.5) and so independent of energy and of mass. Further we assume that a{E, A) 
varies linearly with the wave vector, 

1 



k = — _ ^2^4 ^ (6) 

nc 

and with the form 

a(E, A) ~ 1.02fc - 0.25 . (7) 

Then lo(E,A) were adjusted to ensure that all measured total reaction cross section values 
are matched by using the function form, Eq. (jSj). The resultant optimized values for the 
parameter /q are presented in Fig. ^ The different patterns and colors of the lines show the 
Iq values for different nuclei. The legend of the graph indicates the nuclei. 

While we have used the partial total cross sections from DWA results for neutron scat- 
tering from all the nine nuclei chosen and at all of the energies indicated, from the sets of 
values that result from the fitting process, the two parameters a and e can themselves be 
expressed by the parabolic functions 

a = 1.29 + 0.00250 E - 1.76 x lO""^ E^ , 

e = -1.47 - 0.00234 E + 4.16 x 10"*^ E^ , (8) 

where the target energy E is in MeV. There was no conclusive evidence for a mass variation 
of them. With a and e so fixed, we then adjusted the values of Iq in each case so that actual 
measured neutron total cross-section data were fit using Eq. (0). Numerical values for Iq from 
that process are presented in a table in Ref. jl3]- The values of Iq increase monotonically 
with both mass and energy and that is most evident in Fig. |21 where the optimal values 
lo{E) are presented as different patterned and colored lines. The set for each of the masses 
(from 6 to 238) are given by those that increase in value respectively at 600 MeV. While 
that is obvious for most cases, note that there is some degree of overlap in the values for 
i9^Au and for ^ospb. 

In Figs. Eland EJ we compare the total reaction cross sections generated using the simple 
functional form and tabled values of Iq, and displayed by the red curves, with those ob- 
tained from calculations made using (^-folding optical potentials Blue lines represent the 
predictions obtained from those microscopic optical model calculations. The experimental 
data nj| are depicted by circles. 



The results for scattering from ^Be, ^^C, ^^O, ^^F, ^'^Al, ^^Ca, '^^Cu, ^^Zr and "^Sn are 
displayed in segments (a) through (i) of Fig. [HI respectively. In segment (a), the data are 
well reproduced by both the (^-folding predictions resulting from the folding with the ^Be 
ground state OBDME found with (0+2)^ spectroscopy and by those obtained with the 
simple functional form method. 
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FIG. 1: The parameter values for Iq in the simple functional form for the p-A reaction cross 
sections. 



Calculated p-^^C reaction cross sections are compared with the experimental data in 
segment (b) of Fig. El The reaction cross sections obtained from (yf-folding calculations are 
in good agreement with the experimental data but only in the energy range above 20 MeV. 
On the other hand, the results obtained from the simple functional method are excellent 
for all energies, replicating the data very well even in the lower energy range from 10 MeV. 
There are two data points, at 61 MeV and at 77 MeV, in disagreement with the calculated 
results however. But, as noted previously these data points should be discounted. 



Predictions for p-^^0 and for p-^^F scattering are compared with the data in segments (c) 
and (d) of Fig. El For p-^^0 case, there are many data points at the energies between 20 to 
40 MeV. Predictions from (/-folding calculations while replicating the data very well at and 
above 25 MeV, overestimate at lower energies. That (/-folding result also underestimates the 
datum at 250 MeV; the sole datum above 50 MeV. In contrast, the results obtained from 
the simple functional method are in excellent agreement with the experimental data at all 
energies. For p-^^F, although ^f-folding calculations reproduce the data very well, the simple 
functional form method gives slightly more accurate predictions. 

Total reaction cross section predictions for p-^^Al and for p-^^Ca are compared with the 
experimental data in segments (e) and (f) of Fig. |21 Again while (/-folding calculations 
reproduce the data quite well to 200 MeV, three data points at 180 to 300 MeV are not 
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Neutron energy (MeV) 



FIG. 2: The values of Iq that fit neutron total scattering cross-section data from the nine nuclei 
considered and for energies between 10 and 600 MeV. The curves portray the best fits found by 
taking a function form for Iq{E). 



matched. The (yf-folding results underestimate them noticeably. But predictions from simple 
functional form replicate the data very well at all energies. One data point at 61 MeV is 
exceptional in the set. With ^°Ca, the folding model approach is not expected to be rehable 
at the energies in the range 10 to 20 MeV, as is the case with ^^C, since for excitation 
energies of that range, both nuclei have clearly discrete spectra. That is true for most light 
mass nuclei but little or no total reaction cross sections have been reported for them. Indeed 
the reaction data from both ^^C and from ^''Ca show rather sharp resonance-like features 
below 20 MeV. Both the ^^-folding calculations and functional form calculations reproduce 
the rest of the ^°Ca data very well. 

For ^^Cu, our predictions at low energies may be slightly too small and the parameter 
sets driven too severely by the sole datum at 30 MeV in the range 20 to 70 MeV. Also the 
data in the range 100 to 300 MeV are quite scattered but the simple functional form gives 
a good average result. 

In segment (h) of Fig. El the predicted total reaction cross sections from p-^°Zr scattering 
are compared with the experimental data. Results from gf-folding calculations are in very 
good agreement with the data although the data value at 30 MeV is overestimated. The 
results obtained from the simple functional form are in excellent agreement with the exper- 
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FIG. 3: Energy dependence of an for proton scattering from (a) ^Be, (b) ^^C, (c) ^^O, (d) ^^F, 
(e) ^^Al, (f) ^'^Ca, (g) ^^Cu, (h) ^^Zr and (i) ^^^Sn. Blue lines represent the results obtained from 
^-folding optical potential calculations while the red curves portray the values obtained by using 
simple functional form. 



imental data at the few energies measured, but the shape is not optimally smooth. Lack of 
data meant that we had to use the (yf-folding values to specify the functional form. That is 
also the case with masses 140, 159, and 181. 

The p-^^'^Sn total reaction cross section results are given in segment (i) of Fig. El where 
the two predictions again are compared with the data. Although not as good as the results 
found for scattering from light mass nuclei, the ^^-folding potential still gives reasonable shape 
prediction. However, the model underestimates the data by 5 to 10%. But predictions from 
the simple functional form model form an excellent reproduction of the data at all energies 
except 61 MeV. This 61 MeV data point is again exceptional being much smaller than other 
data and the predictions as in the cases of ^^C and ^'^Al. 

Predictions for p-^^°Ce, p-^^^Tb, p-^^^Ta, p-^^^, p-^*^^ and for p-^^^U scattering are com- 
pared with the (limited amount of) data in segments (a), (b), (c), (d), (e) and (f) of Fig. HI 
respectively. The (^-folding calculations give very good agreement with that data for p-^^°Ce, 
slightly underestimate the data for p-^^^Tb, and for p-^^^Ta, underestimate data at the ener- 
gies to 20 MeV and overestimate data in the energy range 40 to 60 MeV. In all cases results 
predicted by the simple functional form are excellent reproductions of the experimental data. 
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FIG. 4: Same as Fig. El but from (a) i^°Ce, (b) ^^^Tb, (c) ^^^Ta, (d) ^^Uu, (e) ^ospb, and (f) 

238^ 



For ^^^Au, the (/-folding optical potential calculations are in good agreement with most 
data; the 29 MeV datum grossly underestimated by the calculations. But that data point is 
also at odds with the energy trend of the other data. Save for that 29 MeV value, the simple 
functional form gives even better predictions of the data. The energy variation of the p-^'^^Pb 
reaction cross sections is shown in segment (e) of Fig.|3]where the predictions from (/-folding 
optical potential calculations and from the simple functional form calculations are compared 
with a fairly extensive set of experimental data. In making the ^f-folding potentials, we have 
used Skyrme-Hartree-Fock wave functions which have been shown to be more realistic 
than simple oscillator model ones. Still such (/-folding calculations underestimate the data 
up to 50 MeV. But simple functional form calculations are in excellent agreement with the 
experimental data, save that data values at 30, 61 and 77 MeV again are exceptional. The 
predictions of total reaction cross sections for p-^^^U scattering from the (/-folding optical 
potential calculations and from using the simple functional form are compared with the few 
data in segment (f ) of Fig. 0] Given the lack of data the two results are virtually identical. 

As a final note regarding many of the exceptional data values so defined in the foregoing, 
Menet et al. [14] argue that there may be a systematic error in the studies reported in those 
experiments. 

The total neutron scattering cross sections generated using the function form for partial 
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total cross sections with the tabled values of Iq and the energy function forms of Eq. (jH)) for 
a and e, are shown in Fig. El They are displayed by the red lines that closely match the 
data which are portrayed by circles. The data that was taken from a survey by Abfalterer et 
al. 3| which includes data measured at LANSCE that are supplementary and additional to 
those published earlier by Finlay et al. 3]. For comparison we show results obtained from 
calculations made using (/-folding optical potentials 6]. Blue lines represent the predictions 
obtained from those microscopic optical potential calculations. Clearly for energies 300 MeV 
and higher, those predictions fail. 

The total cross sections for neutrons scattered from the nuclei considered are compared 
with data in Fig. El Again the (/-folding potential results are displayed by the blue curves 
while those of the function form are shown by the red curves. We show in segment (h) of 
Fig. the results for neutron scattering from ^"^Pb. In this case we used Skyrme-Hartree- 
Fock model (SKM*) densities to form the (/-folding optical potentials. That structure 
when used to analyze proton and neutron scattering differential cross sections at 65 and 200 
MeV gave quite excellent results Indeed those analyzes were able to show selectivity for 
that SKM* model of structure and for the neutron skin thickness of 0.17 fm that it proposed. 
Using the SKM* model structure, the (/-folding optical potentials gave the total cross sections 
shown by the blue curve in segment (h) of Fig. El Of all the results, we believe these for ^osp]-, 
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point most strongly to a need to improve on the (yf-folding prescription as is used currently 
when energies are at and above pion threshold. Nonetheless, it does do quite well for lower 
energies, most notably giving a reasonable account of the Ramsauer resonances below 
100 MeV. However, as with the other results, these (^-folding values serve only to define a set 
of partial cross sections from which an initial guess at the parameter values of the function 
form is specified. With adjustment that form produces the red curve shown in segment (h) 
of Fig. El which is an excellent reproduction of the data, as it was designed to do. But 
the key feature is that the optimal fit parameter values still vary smoothly with mass and 
energy. 

Without seeking further functional properties of the parameters, one could proceed as we 
have done this far but by using many more cases of target mass and scattering energies so 
that a parameter tabulation as a data base may be formed with which any required value 
of total scattering cross section might be reasonably predicted (i.e. to within a few percent) 
by suitable interpolation on the data base, and the result used in Eq. (0). 

IV. CONCLUSIONS 

Measured reaction cross sections for 10 to 300 MeV proton scattering from nuclei ranging 
in mass from ^Be to ^^^U are well reproduced by calculations made using the g-folding 
model of the optical potential. Measured total cross sections for 10 to 600 MeV neutron 
scattering from nuclei ranging in mass from ^Li to ^38^ are well reproduced. Those calculated 
results gave a set of partial reaction cross sections in each case that vary smoothly with 
target mass and energy. Those variations are well reproduced by a simple three parameter 
function. Using the simple function form total and total reaction cross sections are well 
reproduced. This simple functional method can be used to estimate cross sections for many 
useful applications, such as, medical radiotherapy. 
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